Frustration induced plateaux in S > 1/2 Heisenberg spin ladders 
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We study the T — magnetization of frustrated two-leg spin ladders with arbitrary value of the 
spin 5*. In the strong rung limit, we use degenerate perturbation theory to prove that frustration 
leads to magnetization plateaux at fractional values of the magnetization for all spins S, and to 
determine the critical ratios of parallel to diagonal inter-rung couplings for the appearance of these 
plateaux. These ratios depend both on the plateau and on the spin. To confirm these results, and to 
investigate the properties of these ladders away from the strong coupling limit, we have performed 
extensive density matrix renormalization group (DMRG) calculations for S < 2. For large enough 
inter-rung couplings, all plateaux simply disappear, leading to a magnetization curve typical of 
integer-spin chains in a magnetic field. The intermediate region turns out to be surprisingly rich 
however, with, upon increasing the inter-rung couplings, the development of magnetization jumps 
and, in some cases, the appearance of one or more phase transitions inside a given plateau. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Mg, 75.30.Kz 



I. INTRODUCTION 



Motivation 



Quantum magnets in high magnetic fields possess a 
rich variety of physical phenomena ranging from the ex- 
istence of fractional magnetization plateau 1,2 or the Bosc- 
Einstcin condensation of magnons 3 to the possible exis- 
tence of the spin-equivalent of a supersolid phase. 4 > 5 > 6 ' 7 > 8 
Of particular interest in this context are spin ladder 
systems. 9 From the theoretical point of view, they con- 
stitute an interesting and non trivial step from ID to 2D. 
Several materials are known to realize spin-ladders. Most 
of these systems have spin S = 1/2, but recently systems 
which arc modeled by higher spins have been found. An 
example for this is BIP-TENO, which is considered to be 
a frustrated 5=1 spin ladder. 10 

In this paper, we consider general two-leg ladders with 
arbitrary values of the spin S. We pay particular at- 
tention to the behavior in the strong-rung limit which 
is amenable to an effective description in terms of per- 
turbation theory Starting from this limit, we investigate 
the effect of frustrating interactions, in particular the ap- 
pearance of additional fractional magnetization plateaux, 
which have already been shown to exist in several frus- 
trated quasi-lD systems such as spin-1/2 ladders, 11,12 
spin-1/2 tubes, 13 spin-1/2 chains with nearest and next- 
nearest neighbor exchange 14 and their generalization to 
arbitrary values of S 15 as well as some models of spin-1 
ladders. 16,20 In a combined analysis using perturbative 
methods and the density matrix renormalization group 
(DMRG), we make exact predictions for the existence, 
the position and the sizes of these frustration induced 
plateaux in the strong rung limit. We then extend the 
analysis beyond the strong rung limit and identify an 
intermediate regime in which additional features, i.e., 
jumps in the magnetization curves and phase transitions 
inside the plateaux, are realized and compare to recent 
findings for particular values of the interactions. 
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FIG. 1: (color online) Graphical representation of the ladder 
model, Eq. (1). 



Note that, since we are considering large values of the 
spin, one may wonder whether a semi-classical approach 
would be applicable. For the triangular lattice, the 1/3 
plateau has indeed been shown to be an essentially clas- 
sical phase of the up-up-down type stabilized in a field 
range by thermal fluctuations or by quantum fluctua- 
tions treated at the level of linear spin wave theory, 17 
and this approach has been extended to a number of 
other plateaux. 18 We believe however that such an ap- 
proach cannot account for most of the plateaux reported 
here for two reasons. First of all, the results depend cru- 
cially on the value of the spin (S=l, 3/2 or 2). Besides, 
and more importantly, most plateaux are 'quantum' in 
the terminology of Hida and Affleck. 19 They correspond 
to phases that have no classical counterpart with up and 
down spins, and a semi-classical approximation in terms 
of fluctuations around a classical state is clearly inappli- 
cable. 



B. Model and Methods 



We study ladder systems with frustrating interactions 
between the rungs depicted in Fig. 1 and described by 
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the Hamiltonian 

H = Jj_ Sj,i ■ Sj t 2 — H(Sj i + 5? 2 ) 

i 

+ J\\ / ](Sj t i ■ 5i+i,i + Si^ ■ ^+1,2) 

i 

+ JxX)(^+M-42 + 4i-^+i,a) (!) 

i 

Here, the 5j jQ! are spin operators acting on the sites 
positioned on rung i = 1, ... ,7V and leg a = 1, 2. Jj_, Jy 
and J x denote the rung coupling, the inter-rung parallel 
coupling and the frustrating diagonal interaction, respec- 
tively, as sketched in Fig. 1, and H is the magnetic field. 
In the following, we allow the spins to be of arbitrary 
magnitude 5 and we choose K = 1. In addition, we work 
in units of the energy where J± = 1. 

Throughout the paper, the results will be discussed in 
terms of the magnetization per rung defined by: 

-''rung . 

where N vnng is the number of rungs. This magnetization 
varies between and 25, and for isolated dimers and at 
zero temperature, it is a stepwise function of the mag- 
netic field and takes the integer values 0, 1, ... , 25. 

If Jj_ 3> J\\ + J x , the behavior of the ladder is governed 
by the physics of single dimers and at finite fields magne- 
tization plateaux at 0, 1, . . . , 25* will appear. In the fol- 
lowing, we will refer to these plateaux as integer plateaux. 
However, the inter-rung coupling will induce fluctuations 
between the dimers and in particular the competition be- 
tween J11 and J x will lead to new features. If the interac- 
tions between the rungs are not too strong, these can be 
captured by an effective model which can be derived us- 
ing degenerate perturbation theory. This effective model 
in the present case is an anisotropic 5 = 1/2 XXZ-chain. 
It can be solved exactly using the Bethe ansatz, lead- 
ing to exact predictions on the positions and the sizes of 
possible plateaux. 

The scope of the present paper is two-fold. Start- 
ing from the strong-rung limit, we discuss the effec- 
tive model. In particular, we identify critical val- 
ues of the magnetic field at which in addition to 
the integer plateaux new, frustration induced plateaux 
are created. These predictions are compared to re- 
sults obtained using the density matrix renormalization 
group (DMRG). 21 ' 22 ' 23 ' 24 We consider systems with open 
boundary conditions (OBC) with up to TV = 139 rungs 
and perform, when necessary, 15 sweeps keeping max- 
imally 1200 density matrix eigenstates. Typically, the 
maximum discarded weight is of the order of 10 -10 . We 
extrapolate our finite-size results for the size of the var- 
ious plateaux to the thermodynamic limit and compare 
them to the predictions from the effective model for sys- 
tems with spin up to 5 = 2. In order to estimate the 
error of our results after extrapolation, we compute nu- 
merically the gap for the XXZ-chain and perform the 



same extrapolation of the finite-size data as for the lad- 
der systems. We find that our extrapolated results for 
the chain agree with the Bethe ansatz up to an absolute 
error of the order of 5 • 10~ 4 . Assuming that the error in 
the extrapolation of the data for the ladder system is of a 
similar size, it is thus possible to compare the numerical 
results with a high precision to the results of the effective 
model. 

In the second part of the paper, we leave the strong 
rung limit and consider cases where Jy + J x > Jj_. For 
Jy + J x i» Jj_, the system is described in terms of a sin- 
gle chain with effective spin 5 c h a in = 25. Therefore, the 
behavior at finite fields is governed by the physics of in- 
teger spin chains. Between the strong rung limit and the 
limit of an integer spin chain, we find a crossover region 
in which the plateaux disappear. In this intermediate re- 
gion where the description in terms of the effective model 
is not valid any longer, we find that the magnetization 
curves possess additional interesting features. In partic- 
ular, jumps in the vicinity of the plateaux are observed. 
The size of the plateaux can become non-monotonic when 
increasing Jy + J x , and it is possible to obtain phase tran- 
sitions inside the plateaux as already reported in Ref. 16 
for the special case Jy = J x . 

The paper is organized as follows. In Sec. II A, we 
derive the effective model and formulate its predictions 
for the existence and the size of additional frustration 
induced plateaux. In Sec. II B, we present our DMRG 
results for ladders with 5=1, 5 = 3/2 and 5 = 2. 
For the lowest lying frustration induced plateaux we per- 
form a careful finite size extrapolation and compare with 
the quantitative predictions of the effective model. In 
Sec. Ill, we leave the strong rung limit and consider the 
magnetization behavior for systems with 5 = 1/2 up to 
5 = 3/2 as obtained by the DMRG calculations. In Sec. 
IV, we finally summarize our findings. 



II. FRUSTRATED LADDERS IN THE STRONG 
RUNG LIMIT 

A. Effective model from degenerate perturbation 
theory 

In this section we derive an effective II? Hamiltonian 
in the strong rung limit J± 3> Jy + J x for general values 
of the spin 5. This limit has previously been considered 
in Ref. 25 for the case of non-frustrated ladders (J x = 0). 
Here, we consider both the inter-rung couplings Jy and 
the frustrating interactions J x as perturbations and start 
from isolated dimers which can be treated exactly. This 
problem has been treated in Rcf. 11 for the case 5=1/2 
and in Ref. 10 for 5 = 1. At the critical values of the 
fields at which the levels of the dimers cross (cf. Fig. 2), 
we apply degenerate perturbation theory and derive an 
effective model for each level crossing. (Note that on each 
rung the degeneracy at all level crossings for the present 
Heisenbcrg Hamiltonian can only be two- fold). 
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FIG. 2: (color online) Level crossings of a single rung for 
spins 5=1. The full line denotes the ground state at differ- 
ent values of the magnetic field with the two critical points 
Hci = J± and Hc 2 — 2 J± at which the ground state changes 
from singlet to triplet and from triplet to quintuplet, respec- 
tively. Note that at the critical fields the separation between 
the energy levels is J±, defining the energy scale for which a 
perturbation theoretical treatment in the strong rung limit is 
expected to work. 
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TABLE I: Numerical values of the first few coefficients Dj (S) 
up to S = 5/2 as obtained from Eq. (6). 

In the strong rung limit in the vicinity of the critical 
fields, the physics of the system is described by the two 
states whose energy levels cross, since the energy separa- 
tion with higher energy levels is of the order of J± so that 
their influence can be neglected. Thus, an effective de- 
scription around the level crossings is possible by a model 
which takes into account two degrees of freedom on each 
rung. This suggests to introduce a. S — 1/2 spin chain. 
Working out the details of the degenerate perturbation 
theory, one indeed finds the resulting Hamiltonian to be 
of the form 

i 

(3) 

where the <j x,v,z are the Pauli matrices acting at rung 
position i. Thus, the resulting effective model is a S = 
1/2 XXZ-chain in a magnetic field H c ff. In the following, 



we will derive the parameters J xy , J z and H e g for the 
general case of spin- S ladders. 

1. Mapping at the j th critical magnetic field 

In order to perform the perturbative treatment, we 
rewrite the Hamiltonian as 

H = H + V 

H o = E J±Si,i ■ S l:2 - H C] J2i ( s ti + S i,i) 

i 

y = J\\ E (Xi • + S ia ■ S i+h2 ) (4) 

i 

Here, Hcj denotes the critical field at which the j th 
level crossing takes place. It is given by = jJ± 

with j G {1, . . .25}. At these values of the fields, the 
ground state of the single rung is a superposition of the 
states \j — 1, m = j — 1} and \j, m = j) where the first 
number refers to the total spin on the rung and the sec- 
ond one to its projection in the z direction. Note that 
the ground states of Hq at He arc product states of 
these doubly-degenerate states on the rungs. For two 
adjacent interacting rungs i and i + 1, it is therefore use- 
ful to express the perturbation operator Vt in the basis 
\ a )i ® with a, 6 G {\j - l,m = j - 1) , \j,m = j}}. 

Since the z component of the total spin is conserved un- 
der the action of V, the expectation value of the opera- 
tors Si^k ■ S'i+i,; (k,l = 1,2) can be reduced to two non 
zero matrix elements: (i) the diagonal matrix elements 
(b\i+i Sik^i+i i \ a )i which can be shown to be 

equal to \ (o|^ (b\ i+1 SfS? +1 \a) i \b) i+1 independently of A: 
and I, where S" ■= 5?j + Sf 2 , and (ii) the off-diagonal 
elements, which do not depend on k and I either and are 
given by 

Dj(S) - \ (jJU (j - l,j - 6 \j - l,j - l) l 

(5) 

with O = <S' l + 1 S , i+1 1 + S i 1 S^ +1 v A closed form can be 
obtained by expressing the rung eigenstates in the basis 
of the individual spins of the rung using Clebsch-Gordan 
coefficients: 

Dj(S) = - cg(mi,TO 2 ; j)cg(m 3 ~ l,m 4 ; j - 1) 

mi +m 2 =j 
m3+i»4=3 

x cg(m 3 ,TO 4 ; i)cg(mi - l,m 2 ; j - 1) 

x y/S(S + l) - mi (mi - 1) 

x y/S(S + l)-m 3 (m 3 -l), (6) 

where cg(mi,m 2 ;j) is the Clebsch-Gordan coefficient 
corresponding to the projection of the rung eigenstate 
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TABLE II: Numerical values of the critical ratios of J x /J\\ beyond which frustration induced plateaux appear around fields 
H C] . 



| j, mi +m2) i onto the product state of spins compos- 
ing the rung \S, mi) i 1 \S, m 2 ) i 2 - mi,rri2, m 3 and m 4 G 
{— 5, . . . , S 1 }. The numerical values of the first few coef- 
ficients Dj(S) for S 1 < 5/2 are displayed in Tab. I 

At each He , the mapping to the anisotropic spin i 
chain Eq. (3) is now achieved by introducing the pseudo- 
spin operators <Xj acting on the states \j — 1, m = j ' — 1} 
and \j, m = j) as follows: 



of [J - 1,J - 1) = 

of li,i> = I 



§li-U-i> 



ot u 



li - i,i - 1) = li, j) 



ar|j-l,j-l} = 



(7) 



<r+ \j,3) 







^ li.i) = li - 1, ./ - 1> 

With this convention, the effective model takes on the 
form of an XXZ-chain in terms of these pseudo-spin op- 
erators. The parameters of this effective Hamiltonian can 
be expressed in terms of the original parameters of the 
ladder system as 



J xy ^{S) 



ADj(S) (Jii - Jy 



2 J z j 



(H-H, 



(8) 



Note that J z does not depend on j while the other param- 
eters J xy and H c $ are different at each He ■ This shows 
that the critical behavior, in particular the appearance 
and the size of plateaux, will depend on the level cross- 
ing j. Note also that we expect the region of validity 
of the effective model to become smaller for larger spins 
and higher lying plateaux since the difference between 



the center of the plateaux defined by H^ s 
level crossing point Hc 4 increases with j. 



and the 



2. Position and size of the frustration induced plateaux 

From the expressions (8), the Heiscnbcrg point at 
which J z = J x y>i translates into the condition 

Jx 



J,, 1 ; 8D 3 (S) 



1 



1 



(9) 



As is known from the Bethe ansatz, 26 at this point a 
gap opens for J z > J xv '\ leading to a plateau in the 
magnetization curve of the original ladder system. Since 
these plateaux are only due to the existence of frustrating 
inter-rung interactions, we refer to them in the following 
as frustration induced plateaux. In Tab. II we present 
the critical values of the ratios J x / J\\ for general spin- 
s' ladder systems at which these plateaux appear. Note 
that the Hamiltonian is symmetric under exchange of J x 
and Jii, so that at ratios larger than 1 another critical 
point is obtained at which the plateaux disappear at the 
inverse of the ratios given in Tab. II. 

The size of these frustration induced plateaux can be 
determined by considering the analytic expression for the 
gap obtained from the Bcthc ansatz. For convenience, wc 
rewrite the effective model as 

H = J xy ' j (S) (SfSf +1 + S?S? +1 + pi (S)S Z S Z +1 ) (10) 

where 

J z J\\ + Jx 



P 3 (S) 



j*y,j(S) 8Dj(5)(J|| - Jx) 



(11) 



The size of the gap for p > 1 is then found to be 
(J x y>i(S) = l) 26 



A(p) = sinh($) ^2 



(-iy 



cosh(n$) ' 



cosh($) = p (12) 



5 




FIG. 3: (color online) (a) DMRG results for the magnetization of a S = 1 ladder with N = 41 rungs in the strong rung limit 
as a function of r = Jy/Jx for Jn + J x = 0.1. The solid lines delimiting the plateaux are the results obtained from the Bethe 
ansatz solution of the effective model (3) for the position and the size of the plateaux, (b) Phase diagram as obtained from the 
Bethe ansatz solution of the effective model. For r = Jii / J x ^ 1 the frustration induced plateaux (FP) are separated from the 
integer plateaux (IP) by Luttinger liquid phases (LL), while for r = 1 the plateaux are connected by jumps. 



This function grows exponentially for p » 1, but be- 
comes linear for p — > oo. Thus, from Eq. (12) we obtain 
two approximate expressions for the size of the fractional 
plateau w\ rac : 

<ac = J||+Jx-16^(5)(J||-J x ) 

ifp>l, i.e., J|| » J x (13) 



<ac = 32^(5) (J x - J,,) exp 




if p w 1, i.e., J x /J| ~ r c , (14) 

where we denote by r c the critical value of J x / J\\ at which 
the plateau opens. 

Note that if A(p) > 0, the plateau opens at the value 
of H given by H e g = in both directions when increasing 
and when decreasing H . Thus, the size of the plateau is 
given by twice the size of the gap in the effective model, 
which has already been taken into account in the above 
expressions. In addition, from the effective model using 
the Bethe ansatz it is possible to deduce the range over 
which the magnetization grows until it reaches the next 
integer plateau. Thus, it is possible to obtain expressions 
for the size of the integer plateaux. We find 

<t +1 = \J\\ - Jx I (Dj(S) + D j+1 (S)) (15) 

For the special case J x = J|| the size of the integer 
plateau is always Jj_, as expected from isolated dimers. 
The same considerations lead to expressions for the crit- 
ical fields delimiting the magnetization curves between 
two integer plateaux. We obtain 

# Clow = J±-AD j (S)\J ll -J x \ (16) 
H Chlgh = 2S (Jj_ - Jy - J x ) - AD 3 {S) | J,, - J x |(17) 



Finally, we would like to mention that the prediction for 
both the integer as well as the fractional plateaux is in 
full agreement with general considerations for quantum 
many-body systems which state that a finite excitation 
gap (and therefore a plateau) is only possible if n(S — 
m) = integer, where n is the periodicity of the ground 
state and m the average magnetization per spin. 27,28 In 
the following section we will compare these predictions 
from the effective model to the results of our DMRG 
calculations. 



B. Frustration induced plateaux for S = 1, 5 = 3/2 

and S = 2 

In this section we present our DMRG results for lad- 
der systems with 1 < S < 2 as a function of J x /J|| 
keeping J x + J\\ = 0.1, which we expect to be suffi- 
ciently small for the systems to be in the strong rung 
limit. In Fig. 3(a) we show our results for a finite sys- 
tem with N — 41 rungs for 5=1. For small values 
of r = J x / J|| , the behavior is characterized by a very 
large plateau at M = 1 and a Luttinger liquid (LL) be- 
fore and after this plateau. For r ks 1, however, two new 
plateaux appear, one at M = 1/2 and the other one at 
M = 3/2. These are the frustration induced plateaux 
predicted by the effective model. For a comparison of 
the predictions in the thermodynamic limit for the size 
and position of these additional plateaux with the finite- 
size DMRG results we have included the boundaries of 
the plateaux as obtained from Bethe ansatz. Although 
the system is rather small, the comparison is excellent. 
However, due to the exponentially slow opening of the 
plateaux at the critical fields it is very difficult to esti- 
mate the critical points in this comparison. We will come 
back later to this point and present a detailed compar- 
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FIG. 4: (color online) (a) Size of some frustration induced 
plateaux for different values of S. The lines correspond to the 
Bethe ansatz result for the effective model and the data points 
are the results of finite size extrapolations of DMRG results. 
We estimate the error of the extrapolation to be of the order of 
the symbol size or smaller, see text, (b) Examples of our finite 
size extrapolation for some selected values of the parameters. 
As can be seen, a linear extrapolation is applicable for these 
systems. 



ison after finite size scaling of the numerical results. In 
Fig. 3(b) we present the magnetic phase diagram as ob- 
tained from Bethe ansatz. Two remarks are in order. 
First, the frustration induced plateaux are of different 
size and open at different critical values of r = J x /J»; 
the plateau at M = 3/2 appears at a significantly smaller 
value of r then the plateau at M = 1/2. Indeed, in order 
to observe the lower frustration induced plateau, rather 
large values of Jy./J\\ are needed, so that we expect that 
only strongly frustrated magnetic ladder compounds will 
possess this plateau. Second, at r = 1 the size of the 
plateau as a function of r has a kink. This is due to the 
symmetry of the system under exchange of J x and Jii . 

Now, we turn to the finite-size extrapolation of the 
numerical data. In Fig. 4(a), we present the results of 
this extrapolation for the size of various frustration in- 
duced plateaux as a function of r for systems with S = 1, 



5 = 3/2, and 5 = 2, and in Fig 4(b) an example of the 
finite size scaling for the M = 1/2 plateau in S = 1 sys- 
tems is given. Due to the symmetry under exchange of J x 
and J|| , we only discuss results for r < 1 in the following. 
As can be seen, in all cases the plateau opens exponen- 
tially and possesses the same size Wf rac = (Jx + J\\) at 
r = 1, independent of the value of 5 in accordance to the 
prediction of the effective model, Eq. (13). The overall 
agreement with the Bethe ansatz results is very good, in 
particular for r approaching 1. However, due to the ex- 
ponentially slow opening of the plateau, it is very difficult 
to identify numerically with a high precision the values of 
the critical fields, but up to this uncertainty, they are in 
good agreement with the predicted values. Note that in 
the parameter region where the curvature is maximal the 
extrapolated results show the largest deviation from the 
Bethe ansatz. This deviation seems to be systematic and 
the obtained size of the plateau is found to be smaller 
than the predicted one. This tendency becomes stronger 
when increasing the value of J x + Ju . 

Additional calculations for some selected values of the 
parameters for 5 = 3/2 and 5 = 2 give further support 
for this picture. Since we do not expect any significant 
deviations from the predictions of the effective model, 
we refrain from presenting a detailed analysis for these 
cases. We therefore conclude this section by confirming 
the validity of the description in terms of the effective 
model for values of J x + J\\ < 0.1 for 5 < 2 and expect it 
to be valid for general values of 5. In the next section we 
increase the strength of these interactions and describe 
our findings beyond the strong rung limit. 



III. BEYOND THE STRONG RUNG LIMIT 

Increasing the value of J x + Ju will eventually lead 
to a situation where the spacing of the energy levels of 
the single dimcrs J± is not large enough any more to 
consider the inter-rung couplings as perturbations. For 
J x + J\\ 2> J±, the system can be described in terms of 
a single chain with effective spin 5 G ff = 25, 13 leading to 
a continuously growing magnetization curve without any 
plateau. Between the strong rung limit and the limit of 
integer spin chains, however, we expect to find a crossover 
region in which the physics is not a priori clear. This re- 
gion has been studied in Ref. 16 for the fully frustrated 
case J|| = J x , and various plateaux have been identified 
over a wide range of Ju + J x ■ Interestingly, within some 
of these plateaux, a first order phase transition has been 
found to take place. It is therefore interesting to investi- 
gate the behavior of the integer and frustration induced 
plateaux in the more general case J\\ ^ J x . In order 
to keep a connection to the previous findings, we choose 
values of r so that frustration induced plateaux exist, 
but which are at the same time different enough from 
r = 1 so that new aspects can come into play. Since, to 
our knowledge, this region of the parameter space has not 
yet been investigated for 5=1/2 two- leg ladder systems, 
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FIG. 5: (color online) (a) DMRG results for the magnetization of a S = 1/2 ladder with N = 39 rungs when changing Ju + J x 
and keeping r = J\\/Jx = 0.8 fixed, (b) Phase diagram as obtained from the DMRG calculations for systems with N — 39 
rungs. The notation [i : j] refers to the value of the total spin of two consecutive rungs i and j as described in the text. Thick 
solid lines indicate the position of jumps in the magnetization curves. 



we will start our analysis with this case. In the following, 
we extend it to higher spins up to S — 3/2 in order to 
capture the changes in the magnetization behavior when 
increasing the value of S. Due to the complexity of the 
calculations for these cases, and since the details of the 
transitions are not a main focus of the present work, we 
refrain from performing an elaborate finite size scaling 
analysis at this point and expect that the main features 
are well captured for systems of the sizes presented. A 
more detailed analysis of the nature of the phase transi- 
tions and of the possible significance of finite-size effects 
for S > 1/2 ladders beyond the strong rung limit is left 
for future studies. 

In Fig. 5 we present our DMRG results for S = 1/2 
ladders with TV = 39 rungs as a function of J x + Jii while 
keeping J x /J\\ = 0.8 fixed. There are several features we 
would like to discuss. Even for the rather large values of 
the inter-rung couplings, the fractional plateau at M = 
1/2 exists and is of a size comparable to the one in the 
strong rung limit. At J x + J\\ ~ 1.5, however, the plateau 
vanishes and the magnetization curve resembles the one 
of a S — 1 spin chain. Between J x + Ju w 1.2 and 
the point at which the plateau vanishes, a jump in the 
magnetization curve is visible which starts at values M « 
1 and comes down when increasing the value of J* + Ju. 
It is remarkable that this jump and the plateau both 
disappear around the same values of the interactions. In 
addition to this jump above the fractional plateau, two 
other jumps appear. One reaching the M = h plateau 
from below and the other going down to M = 0. Such 
discontinuities indicate the position of a first order phase 
transition which has been discussed in Ref. 29 to connect 
a region populated by a mixture of singlets and triplets 
with a region consisting only of triplets in the case of the 
jump above the plateau. We expect the other jumps to 
be of similar nature and to delimit the region in which 
the system gets effectively described in terms of a spin-1 



chain against the region in which the ladder physics is 
predominant. The plateau at M = 1/2 is characterized 
by an alternation between rungs in a singlet and a triplet 
state. To describe this, we adapt the notation introduced 
in Ref. 16 and denote possible alternating order by [i : j], 
where the integers i and j denote the values of <S^ ot j = 
Si.i + £j,2 on two consecutive rungs i and j. Note that 
for r = 1, the rungs arc in exact eigenstates of (Sl otal ) 2 , 
while for r / 1 this, in general, is only approximately 
true. In our notation we then choose the value of S£ otal 
which is closest to the next integer value. 

In Fig. 6 we present the magnetic behavior of S = 1 
ladders beyond the strong rung limit. The overall im- 
pression is similar to the one obtained for S — 1/2 sys- 
tems. In particular, magnetization jumps are obtained 
in the vicinity of the plateaux and for J x + Ju > 1.3 
the magnetization curve resembles that of a spin-2 chain. 
The fractional plateaux vanish slightly before the inte- 
ger plateau. For the fractional plateau at M = 3/2, 
and also for the integer plateau at M = 1, the point at 
which the plateaux disappear seems to be connected to 
the existence of the jumps in the magnetization curve. 
An additional interesting feature is obtained for the inte- 
ger plateau at M = 1. At J x + Ju « 1, there is a kink in 
the size of the plateau. At this point, we find indications 
for a phase transition taking place from a phase in which 
all rungs are in the triplet state to a phase with broken 
translational symmetry in which singlets alternate with 
quintuplets. This is shown in more detail in Fig. 6(c), 
where the local magnetizations on the two rungs at the 
center of the system are shown. The plateau at M = 1 
has a uniform spin density up to the critical value of 
J x + J\\ ~ 1 where an alternating pattern between sin- 
glets and quintuplets forms and then suddenly disappears 
when the plateau closes. This transition corresponds to 
the first-order transition found in Ref. 16 for the special 
case J x = J\\ , where at the transition the total spin of the 



8 




J x + J 



FIG. 6: (color online) (a) DMRG results for the magnetization 
of a S = 1 ladder with N = 39 rungs as a function of Ju + J x 
for r = J\\/Jx = 0.9. (b) Phase diagram as obtained from 
the DMRG results for this finite system. The notation [i : j] 
refers to the value of the total spin of two consecutive rungs 
i and j as described in the text. Thick solid lines indicate 
the position of jumps in the magnetization curves, (c) Value 
of (S z ) for the sites of two consecutive rungs at the center of 
the system as a function of Ju + J x inside the plateaux at 
M = 1/2 (O), M = 1 (A) and M = 3/2 (□). 



rungs jumps abruptly from 1 to resp. 2 on alternating 
rungs. In the present case, despite an attempt at a finite- 
size analysis, we have not been able to decide whether the 
transition remains discontinuous, or whether it turns into 
a continuous phase transition of the Ising type, a plau- 
sible alternative in view of the dimerized nature of the 
[2 : 0] phase, and the nature of this transition as well 
as of the disappearance of the plateaux is left for future 



investigation. Note, however, that such a transition from 
a uniform to a symmetry broken state can be expected 
due to the fact that upon increasing the inter-rung cou- 
pling a state in which adjacent spins become more and 
more different is favored since it minimizes the energy 
of the inter-rung couplings. For completeness, we have 
also shown in Fig. 6(c) the magnetization pattern of the 
fractional plateaux. As expected from the analysis in the 
strong rung limit, the fractional plateaux at M = 1/2 
and 3/2 have an alternating pattern between singlet and 
triplet rungs, or between triplet and quintuplet rungs, 
respectively. 

In Fig. 7 finally we present our results for the mag- 
netization curves for the 5 = 3/2 ladder systems. The 
situation here is much richer than for the ladders with 
5 < 1. Again, for values of the inter-rung couplings 
large enough (J x + Jy > 1.25) the magnetization curve 
resembles that of an integer spin chain with 5 = 3 in this 
case. Jumps are visible above the plateaux at M > 1. 
The size of these plateaux becomes non-monotonic when 
increasing the values of the inter-rung interactions, indi- 
cating phase transitions similar to the one found in the 
M = 1 plateau for the 5 = 1 ladder. Note that in the 
region 1.15 < J x + J\\ < 1-25 the lower parts of the mag- 
netization curves show anomalous behavior, like a kink 
around M « 0.6. Despite significant efforts to improve 
the convergence of the DMRG, in this region it turns out 
that the DMRG tends to get stuck at excited states and 
it is difficult to find the true ground states. We depict the 
corresponding magnetization curves in Fig. 7(a) with a 
dotted line. Repeated computation of these curves with 
much stricter convergence parameters shows that, over 
a wide range, the magnetization curves are well repro- 
duced. However, in the interesting regions the quality of 
the calculations remains unclear and we consider these 
results with caution. 

In Fig. 7(c) we analyze the magnetization pattern in- 
side the various plateaux. As expected from the effective 
model, when the inter-rung coupling is not too large, the 
magnetization patterns inside the fractional plateau at 
M = 1/2, 3/2 and 5/2 show a breaking of the trans- 
lational symmetry and the local magnetizations alter- 
nate between values close to the ones of rung-singlets and 
triplets (M = 1/2), triplets and quintuplets (M = 3/2) 
and between quintuplets and septuplets (M = 5/2). As 
for the spin-1 case, and in agreement with the analy- 
sis of Ref. 16, the integer plateaux undergo a transition 
from [1 : 1] to [2 : 0] and from [2 : 2] to [1:3] respectively. 
However, this is not the whole story, and Fig. 7(c) reveals 
additional and unexpected phase transitions. First of all, 
phase transitions are not only found in integer plateaux, 
but there is clear evidence of transitions in the 1/2 and 
3/2 plateau. In the case of the 3/2 plateau, it is quite sim- 
ilar to the tendency observed in Ref. 16, with positive but 
more strongly alternating magnetizations upon increas- 
ing the inter-rung coupling. By contrast, the transition 
that takes place in the 1/2 plateau is to a " fcrrimagnetic" 
state where the magnetization alternates between posi- 
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tive and negative values. Finally, the integer plateau at 
M = 1 undergoes two transitions. After the expected 
transition from the [1 : 1] phase to the [2 : 0] phase, it 
undergoes a second transition to another " ferrimagnetic" 
phase with alternating positive and negative magnetiza- 
tion. 

In order to provide further support for this picture, 
we show in Fig. 8(a) the magnetization on two consecu- 
tive rungs for different system sizes up to N = 99 rungs 
for the M = 1 plateau. Additional evidence is obtained 
by considering the value of the total spin Sl otal on two 
consecutive rungs in the bulk which is done by comput- 
ing the expectation value ((^otai) 2 )- When r = 1, we 
find the rungs to be in exact eigenstates of (S'totai) 2 f° r 
the various system sizes under consideration. In the case 
r = 0.9, we perform a finite size scaling by considering 
system sizes ranging from N = 39 up to N = 99 rungs. 
The results are shown in Fig. 8(b). 

In the case r = 1, we clearly identify three phases real- 
ized on the M = 1 plateau in which the two consecutive 
rungs are in eigenstates of (^otai) 2 with values of the to- 
tal spin [1 : 1], [0 : 2], and [1 : 3]. The first transition 
happens exactly at Ju + J x = 1 as shown in Ref. 16. 
It is a fist order transition connecting a state without 
broken translational symmetry with a state with broken 
translational symmetry The second transition point is 
obtained using numerical results for the ground state en- 
ergies in the third phase, which is described by a chain 
of alternating spin-1 and spin-3 sites and whose ground 
state energy is a linear function of J\\ + J x • The intersec- 
tion point with the ground state energy of the system in 
the second phase, which is a constant since we are deal- 
ing with a product wave function of rung singlets and 
rung quintuplets, results in the critical point which we 
locate at Ju + J x « 1.17. This is fully consistent with 
the numerical results for Sf and Sl otai . 

For the case r = 0.9, after extrapolating to the ther- 
modynamic limit, our results for Sl otal provide essentially 
the same picture. However, the rungs are not in exact 
eigenstates of (•S'totai) 2 ' leading to non-integer effective 
values of the computed 5^ ota j. The data presented in 
Fig. 8(b) indicates the existence of a jump in Sl oted at 
both transitions, supporting the picture that both could 
be of first order. 

The phase transitions to the "ferrimagnetic" configu- 
rations raise an interesting issue. First of all, we note 
that no translational symmetry is broken at the transi- 
tion since all phases are dimerized. But at the same time, 
since the total spins of the rungs are not conserved when 
r^l, these operators cannot be used as conserved quan- 
tities to characterize the phases. It thus remains an open 
problem to find a way to characterize these phases. This 
is reminiscent of the problem of characterizing different 
singlet phases in spin-1 chains with additional interac- 
tions, or in spin-1/2 ladders. In that context, non-local 
string order parameters have been shown to have different 
values in different phases, and to provide the appropriate 
way to distinguish the phases. 30,31 ' 32 Whether such non- 
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FIG. 7: (color online) (a) DMRG results for the magnetization 
of a 5* = 3/2 ladder with N — 39 rungs as a function of Jii + J x 
for r = J\\/Jx = 0.9. (b) Phase diagram as obtained from the 
DMRG results for these finite systems. The notation [i : j] 
denotes the value of the total spin of two consecutive rungs 
i and j as described in the text. Thick solid lines indicate 
the position of jumps in the magnetization curves, (c) Value 
of (S z ) on the sites on two consecutive rungs at the center 
of the system as a function of J« + J x inside the plateaux at 
M = 1/2 (O), M = 1 (x), M = 3/2 (□), M = 2 (*) and 
M = 5/2 (A). 



local order parameters can be constructed in the present 
case is left for future investigation. 
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IV. SUMMARY AND CONCLUSION 

In conclusion, we have investigated the magnetic be- 
havior of general Heisenberg spin-5 two-lcg-laddcrs in a 
magnetic field. Starting from the strong rung-limit, we 
describe the physics of the systems (independent of the 
value of 5) in terms of an effective 5 — 1/2 XXZ-chain 
obtained from degenerate perturbation theory for values 
of the interactions J x +J\\ <C J±. In this limit, we predict 
the existence of additional fractional plateaux which are 
purely frustration induced and provide exact values for 
the position and the size of these plateaux. Within the 
accuracy of our numerical resolution, we confirm these 
predictions with our DMRG calculations and find that 
the effective model for 5 < 2 is a qualitatively good 
description up to values of J x + Jn rs J±. In the op- 
posite limit, for J x + Jn 3> J±, the magnetizations are 
reminiscent of the ones of integer 25 spin chains. In an 
intermediate regime around J x + Jn = 1.1 Ji to 1.5Jj_ 
(depending on the actual value of 5) , we find additional 
interesting features in the magnetization curves. Of par- 
ticular significance are jumps in the vicinity of some of 
the plateaux which can also be realized at M = and 



phase transitions inside the plateaux which can be visi- 
ble as kinks in the size of the plateaux. Some of these 
phase transitions were predicted in Ref. 16 for the spe- 
cial case J x = J\\ for general values of 5 and were found 
to be of first order. In this special case, plateaux at 
M = 0, 1/2 and at M = 25, 25-1/2 arc found not to 
possess these phase transitions, while the others should 
all possess at least one first order transition inside the 
plateaux. Although we leave the question on the nature 
of the transitions open, our findings for finite systems 
for 5 < 1 are in qualitative agreement with the r = 1 
case. However, we find indications for additional, possi- 
bly first-order transitions inside the M = 1 and M = 1/2 
plateaux in the 5 = 3/2 case into 'ferrimagnetic' phases 
where the rung magnetizations alternate between posi- 
tive and negative values. At these transitions, no sym- 
metry is broken. Whether non-local order parameters 
can be devised to characterize these phases is left for fu- 
ture investigation. 
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